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Abstract. IceCube, a future km 3 antarctic ice Cherenkov neutrino telescope, is 
highly sensitive to a galactic supernova (SN) neutrino burst. The Cherenkov light 
corresponding to the total energy deposited by the SN neutrinos in the ice can be 
measured relative to background fluctuations with a statistical precision much better 
than 1%. If the SN is viewed through the Earth, the matter effect on neutrino 
oscillations can change the signal by more than 5%, depending on the flavor-dependent 
source spectra and the neutrino mixing parameters. Therefore, IceCube together with 
another high-statistics experiment like Hyper-Kamiokande can detect the Earth effect, 
an observation that would identify specific neutrino mixing scenarios that are difficult 
to pin down with long-baseline experiments. In particular, the normal mass hierarchy 
can be clearly detected if the third mixing angle is not too small, sin 2 #13 > 10 -3 . 
The small flavor-dependent differences of the SN neutrino fluxes and spectra that are 
found in state-of-the-art simulations suffice for this purpose. Although the absolute 
calibration uncertainty at IceCube may exceed 5%, the Earth effect would typically 
vary by a large amount over the duration of the SN signal, obviating the need for 
a precise calibration. Therefore, IceCube with its unique geographic location and 
expected longevity can play a decisive role as a "co-detector" to measure SN neutrino 
oscillations. It is also a powerful stand-alone SN detector that can verify the delayed- 
explosion scenario. 



1. Introduction 

The antarctic neutrino telescope AMANDA ^ |2| and the future km 3 IceCube E] 
are designed to observe high-energy neutrinos from astrophysical sources. The ice is 
instrumented with photomultipliers to pick up the Cherenkov light from secondary 
charged particles. In order to reach the large volume needed to detect the expected 
small fluxes at high energies, the density of optical modules is far too sparse to measure, 
for example, solar neutrinos. However, it has been recognized for a long time that these 
instruments can detect a supernova (SN) neutrino burst because the Cherenkov glow 
of the ice can be identified as time-correlated noise among all phototubes |H| • This 
approach has been used by AMANDA to exclude the occurrence of a galactic SN over 
a recent observation period [Zj. 

For AMANDA the physics potential of a possible SN observation is essentially 
limited to its detection, notably in the context of the Supernova Early Warning System 



Supernova Neutrino Oscillations at IceCube 



2 



(SNEWS) that would alert the astronomical community several hours before the optical 
explosion jHJEl- F° r the future IceCube with 4800 optical modules, however, the number 
of detected Cherenkov photons would be of order 10 6 and thus so large that several 
interesting physics questions could be addressed in earnest. 

The observed quantity is the number of Cherenkov photons caused by the SN 
neutrinos as a function of time, i.e. a measure of the energy deposited by the neutrinos 
in the ice. Therefore, the information about the SN signal is far more limited than 
what can be extracted from a high-statistics observation in Super-Kamiokande or other 
low-energy experiments that detect individual events. However, galactic SNe are so 
rare, perhaps a few per century, that the chances of observing one depend crucially on 
the long-term stability of the neutrino observatories. Once IceCube has been built it 
may well operate for several decades, backing up the low-energy experiments. Besides 
the detection and associated early warning one could measure important details of the 
neutrino light-curve, for example the existence and duration of the initial SN accretion 
phase, the overall duration of the cooling phase, and so forth. Such an observation 
would provide a plethora of astrophysically valuable information. 

However, from the perspective of neutrino physics a simultaneous observation 
in both IceCube and another large detector such as Super-Kamiokande or Hyper- 
Kamiokande would be especially useful. Assuming that the neutrinos have traversed 
significantly different paths through the Earth, the two signals could well show 
measurable differences caused by neutrino oscillations in matter [TUJ HU H2 ECU El- 
As this Earth effect shows up only for certain combinations of neutrino mixing 
parameters, a dual observation may well distinguish, for example, between the normal 
and inverted neutrino mass hierarchy. It is well known that observing SN neutrinos 
with two or more detectors with different Earth-crossing lengths is extremely useful, 
but IceCube's potential has not been explored in this context. With all the low-energy 
observatories being in the northern hemisphere, IceCube's location in Antarctica is 
uniquely complementary for this purpose. 

Any oscillation signature depends on the small flavor-dependent differences between 
the fluxes and spectra at the source. If these differences were as large as had been 
assumed until recently there would be little question about IceCube's usefulness for co- 
detecting the Earth effect. However, a more systematic study of the flavor-dependence 
of the SN neutrino fluxes and spectra reveals that these differences are more subtle, 
although by no means negligible [13 EH H3 EE] ■ We evaluate IceCube's potential as a 
co-detector from the perspective of these "pessimistic" assumptions about the primary 
fluxes and spectra. 

This paper is organized as follows. In Sec. we show that the neutrino signal 
from a galactic supernova can be measured at IceCube with a sub-percent statistical 
precision. In Sec. 01 we calculate the Earth matter effects on this signal and illustrate 
that it is possible to detect them in conjunction with another high statistics experiment. 
Sec. |U concludes. 
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2. Supernova Neutrino Detection in Ice Cherenkov Detectors 

2.1. Cherenkov Photons in One Optical Module 

The SN neutrinos streaming through the antarctic ice interact according to v e p — > ne + 
and some other less important reactions. The positrons, in turn, emit Cherenkov light 
that is picked up by the optical modules (OMs) frozen into the ice. While the expected 
number of detected photons per OM was calculated in Refs. E], we revisit their 
estimate for two reasons. First, the SN signal was directly scaled to the historical 
SN 1987A observation in Kamiokande II so that the exact assumptions about the 
neutrino flux are not directly apparent. Second, the expected number of Cherenkov 
photons detected by one OM was based on estimating an effective ice volume seen by 
one OM. However, it is much simpler to work in the opposite direction and start with the 
homogeneous and isotropic Cherenkov glow of the ice caused by the SN neutrinos. The 
OM is immersed in this diffuse bath of photons and picks up a number corresponding 
to its angular acceptance and quantum efficiency 

As a first simplification we limit ourselves to the signal caused by the inverse /3 
reaction v e p — > ne + . The v e fluence (time- integrated flux) at Earth is 

T Pe = 1.745 xlO 11 cm" 2 / SN . (1) 

We define the "SN fudge factor" as 

, £ ge ,tot 15MeV / I0kpc \ 2 

/SN " 5 x 10 52 erg (E Pe ) { D ) ' U 

where E Peyto t is the total energy leaving the SN in the form of v e after flavor oscillations 
have been included, (E Pe ) is the average v e energy, and D the distance. 

The energy deposited in the ice per target proton is T Pe (E Pe a). For the inverse f3 
cross section we ignore weak-magnetism and recoil corrections and also the difference 
between v e and positron energy so that jTH] 

, = 9.52X10- cm' (^)\ (3) 

For the neutrino flux of each neutrino and anti-neutrino species we assume a distribution 
of the form 



$ (l + a) 1+a ( E 



(oc + 1) — 

£/ 



(4) 



where E is the average energy, a a parameter that typically takes on values 2.5-5 
depending on the flavor and the phase of neutrino emission, and $o the overall flux at 
the detector in units of cm -2 s _1 . This distribution implies 

{ED = (3 | 1 ^ ) t a) W = £ (E,,f for a = 3 . (5) 

Altogether we thus find 

(E Pe a) = 6.024 x 10" 40 MeV cm 2 f a (6) 
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with 

_ 8(3 + a)(2 + g) / (E Pe ) V 
Ja ~ 15(1 + a) 2 ^15 MeV/ ' {) 

This fudge factor can also be taken to include deviations from the simplified energy 
dependence of the cross section and deviations from the assumed spectral shape. 

The Cherenkov angle for photon emission by a charged particle is cos = (n/3) _1 
where n is the medium's refractive index and (3 the particle's velocity. With n = 1.31 
for ice, neglecting the A-dependence, and j3 — 1 we have = 40.2°. A particle with 
unit charge produces Cherenkov photons per unit path length and per unit wavelength 
band according to 

rf 2 iV 7 _ 2vra sin 2 

dx~d\ ~ A 2 " ' ( ' 

where a = 1/137 is the fine-structure constant. Assuming that n and thus are 
independent of wavelength we integrate over A and find 

d N -y °° „„„ 1 300 nm 

— ^ = 638 cm" 1 . 9 

ax \ A 

Taking the useful wavelength range to be 300-600 nm this translates into 319 photons 

per cm pathlength. Taking the positron mean free path to be 12 cm for an energy of 

20 MeV, and taking it to be proportional to its energy, the number of useful Cherenkov 

photons per deposited neutrino energy is 

^ = 191 MeV" 1 /ch (10) 

with yet another fudge factor fch- 

The density of ice is 0.924 g cm -3 , corresponding to about 6.18 x 10 22 cm -3 
proton targets. Therefore, the SN neutrinos produce 1.241 x 10~ 3 cm -3 / SN f a f C h 
useful Cherenkov photons per unit volume of ice. Multiplying this number with the 
speed of light and dividing by 4tt gives us the resulting diffuse photon flux in units 
of cm -2 s -1 ster -1 . However, the average lifetime of these photons is ci? a b s with R^hs 
the absorption length. Therefore, the neutrino-induced photon fluence is found by 
multiplying the flux with c.R a bs, 

= 0.9874 cm" 2 ster" 1 / SN /„ /ch /abs (11) 

where / abs = i? abs /100m. 

The number of events produced by this fluence in a given OM depends on the 
average quantum efficiency taken to be Q = 0.20. In addition, it depends on the 
angular acceptance, i.e. the effective photo cathode detection area A cat times the angular 
acceptance range f2 acc - Therefore, in one OM we expect 

^events = 310 /sn fa fch /abs foM (12) 

with 

Q A ca t ^acc fT Q\ 

_ bT20 2^b _ cm 2 ~2n~ ' [U} 
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This result is independent of the presence of bubbles in the ice that scatter the photons. 
The Cherenkov glow of the ice represents an isotropic and homogeneous distribution 
that is not changed by elastic scattering. 

2.2. Comparing With Previous Work 

In order to compare our result with the one derived in Ref. |H| we need to translate their 
assumptions into our fudge factors. The v e distribution was taken to follow a Fermi- 
Dirac spectrum with T = 4 MeV, implying (E Pe ) = 12.61 MeV. The distance of the SN 
was taken to be 10 kpc, and the total energy release was scaled to the Kamiokande II 
signal for SN 1987A. With our choice of the j3 cross section these assumptions correspond 
to Ep e>tot = 3.17 x 10 52 erg, i.e. to /sn = 0.754. These authors also used a quadratic 
energy dependence of the cross section. Integrating over their Fermi-Dirac spectrum 
they effectively used f a = 0.663. Further, they assumed 3000 useful Cherenkov photons 
for 20 MeV deposited energy, i.e. effectively fch = 0.785. For the absorption length they 
used 300 m, i.e. / a b s = 3. Finally, they assumed a quantum efficiency of 25%, a cathode 
area of 280 cm 2 , and an acceptance range of 2n, i.e. /om = 1.12. Altogether, we find for 
these assumptions N events = 409 per OM. This compares with 273 in Ref. 0, i.e. our 
result is larger by a factor 1.5. 

The result in Ref. jnj was backed up by a detailed Monte Carlo treatment of 
the production and propagation of Cherenkov photons in the AMANDA detector. 
Therefore, the difference may well relate to details of the OM acceptance and 
wavelength-dependent quantum efficiency and photon propagation. Many of these 
details will be different in IceCube where 10 inch photomultiplier tubes and different 
regions of ice will be used. Detailed values for the detector-dependent fudge factors 
must be determined specifically for IceCube once it has been built. The main difference 
between the assumptions in Ref. jHj and our estimate is the absorption length. When 
using AMANDA as a SN observatory a realistic value was taken to be around 100 m [7j. 
The vast difference between these estimates is that the former was based on the measured 
absorption length in a dust-free region of the ice. For our further estimates we stick to 
100 m as a conservative assumption. 

2.3. Event Rate vs. Neutrino Luminosity 

In our derivation we have used the time-integrated neutrino flux, amounting to the 
assumption of a stationary situation. The absorption time for photons is very small, 
Tabs = -Rabs/c = 0.33 /is i? a b s /100 hi. The SN signal will vary on time scales exceeding 
10 ms. Therefore, the Cherenkov glow of the ice follows the time-variation of the SN 
signal without discernible inertia. Hence one may replace the neutrino fluence with a 
time-dependent flux and N events with an event rate r events . 

Moreover, for our further discussion it will be useful to consolidate our fudge-factors 
into one describing the detector response, and others characterizing the neutrino flux. 
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Therefore, we summarize our prediction for the event rate per OM in the form 

Tevents = 62 S 1 — — — r ( ) /flux /dot (14) 

10 52 erg s _1 \ D / 

where L Pf , is the z/ e luminosity after flavor oscillations and 

15 MeV 8 (Ej) 
/flux (£ Pc > 15 (15 MeV) 3 ' K ' 

100 m O20 250 cm 2 Tvr 
Here, /det also includes corrections for the energy dependence of the (3 cross section. 

We stress that our simple estimate of the counting rate primarily serves the purpose 
of determining its magnitude relative to the background. The important feature is that 
the signal relative to the background can be determined with a good statistical precision. 
Of course, for an absolute detector calibration a detailed modeling would be necessary. 
For our present purpose, however, even an uncertainty of several 10% in our estimated 
counting rate is irrelevant. 



■^cat J 'acc / „ \ 



2.4- Supernova Signal in IceCube 

IceCube will have 4800 OMs so that one expects a total event number of 1.50 x 10 6 , 
taking all fudge factors to be unity. Assuming a background counting rate of 300 Hz 
per OM over as much as 10 s this compares with a background rate of 1.44 x 10 7 . 
Assuming Poisson fluctuations, the uncertainty of this number is 3.8 x 10 3 , i.e. 0.25% of 
the SN signal. Therefore, one can determine the SN signal with a statistical sub-percent 
precision, ignoring for now problems of absolute detector calibration. 

In order to illustrate the statistical power of IceCube to observe a SN signal we 
use two different numerical SN simulations. The first was performed by the Livermore 
group [201 ^ na t involves traditional input physics for mu- and tau-neutrino interactions 
and a flux-limited diffusion scheme for treating neutrino transport. The great advantage 
of this simulation is that it covers the full evolution from infall over the explosion to 
the Kelvin-Helmholtz cooling phase of the newly formed neutron star. We show the 
Livermore v e and v x lightcurves in Fig ^ (left panels). Here and in the following we take 
v x to stand for either z/ M or v T . Apart from very small differences the SN fluxes and 
spectra are thought to be equal for z/ M , u T , z/ M , and v T . 

Our second simulation was performed with the Garching code |21] . It includes 
all relevant neutrino interaction rates, including nucleon bremsstrahlung, neutrino 
pair processes, weak magnetism, nucleon recoils, and nuclear correlation effects. The 
neutrino transport part is based on a Boltzmann solver. The neutrino-radiation 
hydrodynamics program allows one to perform spherically symmetric as well as multi- 
dimensional simulations. The progenitor model is a 15 M© star with a 1.28 M iron core. 
The period from shock formation to 468 ms after bounce was evolved in two dimensions. 
The subsequent evolution of the model is simulated in spherical symmetry. At 150 ms 
the explosion sets in, although a small modification of the Boltzmann transport was 
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necessary to allow this to happen [22] • Unmanipulated full-scale models with an accurate 
treatment of the microphysics currently do not obtain explosions [2H] ■ This run will be 
continued beyond the current epoch of 750 ms post bounce; we here use the preliminary 
results currently available [18]. We show the Garching v e and v x lightcurves in Fig^ 
(right panels). 

We take the Livermore simulation to represent traditional predictions for flavor- 
dependent SN neutrino fluxes and spectra that were used in many previous discussions 
of SN neutrino oscillations. The Garching simulation is taken to represent a situation 
when the v x interactions are more systematically included so that the flavor-dependent 
spectra and fluxes are more similar than had been assumed previously [TH| ITT)| ITT] IT%] . 
We think it is useful to juxtapose the IceCube response for both cases. 

Another difference is that in Livermore the accretion phase lasts longer. Since the 
explosion mechanism is not finally settled, it is not obvious which case is more realistic. 
Moreover, there could be differences between different SNe. The overall features are 
certainly comparable between the two simulations. 

In Fig. El we show the expected counting rates in IceCube on the basis of Eq. (JUj) 
for an assumed distance of 10 kpc and 4800 OMs for the Livermore (left) and Garching 
(right) simulations. We also show this signal in 50 ms bins where we have added noise 
from a background of 300 Hz per OM. The baseline is at the average background rate 
so that negative counts correspond to downward background fluctuations. 

One could easily identify the existence and duration of the accretion phase and 
thus test the standard delayed-explosion scenario. One could also measure the overall 
duration of the cooling phase and thus exclude the presence of significant exotic 
energy losses. Therefore, many of the particle-physics limits based on the SN 1987A 
neutrinos [21] could be supported with a statistically serious signal. If the SN core were 
to collapse to a black hole after some time, the sudden turn-off of the neutrino flux could 
be identified. In short, when a galactic SN occurs, IceCube is a powerful stand-alone 
neutrino detector, providing us with a plethora of information that is of fundamental 
astrophysical and particle-physics interest. 

In addition, IceCube is extremely useful as a co-detector with another high-statistics 
observatory to measure neutrino oscillation effects, a topic that we now explore. 
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Figure 1. Supernova v e and v x light curves and average energies. Left: Livermore 
simulation [201 - Right: Garching simulation |18j . 
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Figure 2. Supernova signal in IceCube assuming a distance of 10 kpc, based on the 
Livermore simulation (left) and the Garching one (right), in both cases ignoring flavor 
oscillations. In the bottom panels we have used 50 ms bins and have added noise from 
a background rate of 300 Hz per OM. 
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3. Observing Supernova Neutrinos Through The Earth 

3.1. Earth Matter Effect on SN Neutrino Spectra 

Neutrino oscillations are now firmly established by measurements of solar and atmo- 
spheric neutrinos and the KamLAND and K2K long-baseline experiments |25[ EEl EZ1 
123 EIIl ED]- Evidently the weak interaction eigenstates u e , and v T are non-trivial 
superpositions of three mass eigenstates u\, z/ 2 and u 3 , 



(17) 



( Ue \ 






Vil 1 






\Vr) 




U 



where U is the leptonic mixing matrix that can be written in the canonical form 



/ 1 



U 







e l5 s 13 








(18) 



C13 

C 23 S 23 I I I I -S|„ f,o 

VO -s 23 c 2 J \-e- iS s 13 ci 3 / V 1. 
Here ci 2 = cos #12 and si 2 = sin6 l i 2 etc., and 5 is a phase that can lead to CP- violating 
effects, that are, however, irrelevant for SN neutrinos. 

The mass squared differences relevant for the atmospheric and solar neutrino 
oscillations obey a hierarchy Am^ tm ^> AniQ. This hierarchy, combined with the 
observed smallness of the angle #13 at CHOOZ [31] implies that the atmospheric neutrino 
oscillations essentially decouple from the solar ones and each of these is dominated by 
only one of the mixing angles. The atmospheric neutrino oscillations are controlled by 
# 23 that may well be maximal (45°). The solar case is dominated by #i 2 , that is large 
but not maximal. From a global 3-flavor analysis of all data one finds the 3a ranges for 
the mass differences Am?- = mf — m' 



and mixing angles summarized in Tabled 



Table 1. Neutrino mixing parameters from a global analysis of all experiments 
(3ct ranges) |2~5*] . 



Observation 


Mixing an 


gle 


Am 2 [meV 2 ] 


Sun, KamLAND 


O12 = 2T- 


-42° 


Amfi = 55-190 


Atmosphere, K2K 


023 = 32°- 


-60° 


|Am§ 2 | = 1400-6000 


CHOOZ 


613 < 14° 




Am§ x w Am§ 2 



A SN core is essentially a neutrino blackbody source, but small flavor-dependent 
differences of the fluxes and spectra remain. We denote the fluxes of u e and v x at Earth 
that would be observable in the absence of oscillations by F° and F°, respectively. In 
the presence of oscillations a u e detector actually observes 

*£, (19) 

where p D (E) is the v e survival probability after propagation through the SN mantle and 
perhaps part of the Earth before reaching the detector. 



F s D (E)=p D (E)F°(E) + 



(E) 
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A significant modification of the survival probability due to the propagation through 
the Earth appears only for those combinations of neutrino mixing parameters shown in 
Table El The Earth matter effect depends strongly on two parameters, the sign of Am§ 2 
and the value of |0i3| [TU1 lllj. The "normal hierarchy" corresponds to mi < m 2 < m 3 , 
i.e. Am| 2 > 0, whereas the "inverted hierarchy" corresponds to 7713 < mi < m2, i.e. 
Am 2 2 < 0. Note that the presence or absence of the Earth effect discriminates between 
values of sin 2 # 13 less or greater than 10 -3 , i.e. #13 less or larger than about 1.8°. Thus, 
the Earth effect is sensitive to values of #13 that are much smaller than the current limit. 

Table 2. The Earth effect appears for the indicated flavors in a SN signal. 

13-Mixing Normal Hierarchy Inverted Hierarchy 

sin 2 #13 < v e and v e v e and v e 

sin 2 #13 > 1CT 3 v e v e 



Let us consider those scenarios where the mass hierarchy and the value of #13 are 
such that the Earth effect appears for v e . In such cases the v e survival probability p D (E) 
is given by 

f « cos 2 12 - sin 20® sm(20® - 20 12 ) sin 2 (l2.5 , (20) 

where the energy dependence of all quantities will always be implicit. Here 9% is 
the mixing angle between v e and i><i in Earth matter while Am^ is the mass squared 
difference between the two anti-neutrino mass eigenstates V\ and z/ 2 in units of 10 _5 eV 2 , 
L is the distance traveled through the Earth in units of 1000 km, and E is the neutrino 
energy in MeV. We have assumed a constant matter density inside the Earth, which is 
a good approximation for L < 10.5, i.e. as long as the neutrinos do not pass through 
the core of the Earth. 



3.2. Magnitude of Observable Effect at IceCube 

In order to calculate the extent of the Earth effect for IceCube, we will assume that the 
relevant mixing parameters are Am 2 2 = 6 x 10" 5 eV 2 and sin 2 (20i 2 ) = 0.9. We further 
assume that the source spectra are given by the functional form Eq. (J3J. The values of 
the parameters a and (E) for both the v e and u x spectra are in general time dependent. 

In Fig. El we show the variation of the expected IceCube signal with Earth-crossing 
length L for the two sets of parameters detailed in Table El The first could be 
representative of the accretion phase, the second of the cooling signal. We use the 
two-density approximation for the Earth density profile, where the core has a density 
of 11.5 g cm -3 and a radius of 3500 km, while the density of the Earth mantle was 
taken to be 4.5 g cm -3 . We observe that for short distances, corresponding to near- 
horizontal neutrino trajectories, the signal varies strongly with L. Between about 3,000 
and 10,500 km it reaches an asymptotic value that we call the "asymptotic mantle 
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value." For Case (a), this value corresponds to about 1.5% depletion of the signal, 
whereas for (b) it corresponds to about 6.5% depletion. 




Earth— Crossing Length [km] 

Figure 3. Variation of the expected IceCube signal with neutrino Earth crossing 
length L for the assumed flux and mixing parameters of Table The signal is 
normalized to 1 when no Earth effect is present, i.e. for L = 0. The dashed line 
is for the case representing the accretion phase, the solid line for the cooling phase. 



Table 3. Flux parameters for two representative cases. 



Example 


Phase 


(Eo.) 


(Ee.) 






$0 /<TjO 


Asymptotic 






[MeV] 


[MeV] 








Earth Effect 


(a) 


Accretion 


15 


17 


4 


3 


1.5 


-1.5% 


(b) 


Cooling 


15 


18 


3 


3 


0.8 


-6.5% 



Beyond an Earth-crossing length of ~10,500 km, the neutrinos have to cross the 
Earth core with another large jump in density. The core effects change the asymptotic 
mantle value by ~ 1% as can be seen in Fig. El We neglect the core effects in the 
following analysis, and the "asymptotic value" always refers to the asymptotic mantle 
value. 

For the largest part of the sky the Earth effect either appears with this asymptotic 
value ("neutrinos coming from below"), or it does not appear at all ("neutrinos from 
above"). Therefore, we now focus on the asymptotic value and study how the signal 
modification depends on the assumed flux parameters. In Table |U we show the signal 
modification for (E Be ) = 15 MeV, a Pe = 4.0, and a Px = 3.0 as a function of (E Px ) and 
the flux ratio <&? /$? . In Table Owe show the same with a 9 = a D = 3.0. The results 
are shown in the form of contour plots in Fig. 0] 

Even for mildly different fluxes or spectra the signal modification is several percent, 
by far exceeding the statistical uncertainty of the IceCube signal, although the absolute 
calibration of IceCube may remain uncertain to within several percent. However, the 
signal modification will vary with time during the SN burst. During the early accretion 
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Table 4. Asymptotic IceCube signal modification by the Earth effect. The fixed flux 
parameters are (E Pc ) — 15 MeV, ap e = 4.0, and cta x = 3.0. 



Flux ratio (E Px ) [MeV] 

$g e /$? x 15 16 17 18 19 20 



2 


.0 


1. 


,026 


1 


,014 


1 


,002 





,988 





.974 


0, 


,960 


1 


.9 


1. 


.023 


1 


,011 


0, 


,999 


0, 


.985 





.971 


0, 


,956 


1 


.8 


1. 


,021 


1 


.009 


0. 


,995 





,982 





.967 


0, 


,952 


1 


.7 


1. 


,018 


1 


.005 


0, 


.992 


0, 


,978 





.963 


0, 


,948 


1 


.6 


1. 


,015 


1 


.002 


0. 


.988 


0. 


,974 





.959 


0, 


,944 


1 


.5 


1. 


,012 





.998 


0, 


.984 


0, 


,969 





.954 


0, 


,939 


1 


4 


1. 


,008 





.994 





,980 





.965 





.949 


0, 


,934 


1 


.3 


1. 


,004 





.990 


0, 


,975 


0, 


,960 





.944 


0, 


,928 


1 


.2 


1. 


,000 





.985 





.970 





,954 





.938 


0, 


,922 


1 


4 


0. 


,995 





.980 


0, 


.964 


0, 


.948 





.932 


0, 


,915 


1 


.0 


0. 


,989 





.974 





,957 


0. 


,941 





.925 


0, 


,908 





.9 


0. 


,983 
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Table 5. Same as Table0]with ap e = a Px = 3.0. 



Flux ratio (E Dx ) [MeV] 

/$? 15 16 17 18 19 20 



2 


,0 


1 


.036 


1 


.024 


1, 


.012 


1, 


,000 





.986 


0. 


,972 


1, 


,9 


1 


.033 


1 


.022 


1, 


,010 


0, 


,996 





.983 


0. 


,968 


1, 


.8 


1 


.031 


1 


.019 


1, 


,006 


0, 


.993 





.979 


0. 


,964 


1 


.7 


1 


.028 


1 


.016 


1, 


.003 


0, 


,989 





.975 


0. 


,960 


1 


.6 


1 


.025 


1 


.013 





,999 





.985 





.971 


0, 


,955 


1 


.5 


1 


.022 


1 


.009 


0, 


,995 


0, 


,981 





.966 


0, 


,951 


1 


.4 


1 


.019 


1 


.005 





.991 





.976 





.961 


0, 


,945 


1 


.3 


1 


.015 


1 


.001 


0, 


,986 


0, 


,971 





.955 


0, 


,940 


1 


.2 


1 


.010 





.996 





,981 





.965 





.949 


0, 


,933 


1 


.1 


1 


.006 





.991 


0, 


,975 


0, 


,959 





.943 


0, 


,927 


1 


.0 


1 


.000 





.985 


0. 


,969 


0. 


.952 





.936 


0, 


,919 





.9 





.994 





.978 


0, 


.961 


0, 


.945 





.928 


0, 


,911 





.8 





.986 





.970 





.953 





.936 





.919 


0, 


,903 





.7 





.978 





.961 


0, 


.944 


0, 


.926 





.910 


0, 


,893 





.6 





.968 





.950 





.933 





,916 





.899 


0, 


,882 





.5 





.956 





.938 


0, 


.920 


0, 


.903 





.886 


0, 


,870 



Supernova Neutrino Oscillations at IceCube 



13 



1.9 

1.7. 
1.5' 
13 
1.1 
0.9 
0.7 



(a) 



15 16 



17 



IceCube Signal 




_ a; 



19 20 



15 



(b) 




16 



17 



18 19 20 



< E 



Figure 4. Asymptotic IceCube signal modification by the Earth effect. The fixed flux 
parameters are (a) (i?p c ) = 15 MeV, a Pc = 4.0, and = 3.0 and (b) (E^J = 15 MeV, 
Q^e = a Px — 3.0. The contours are equally spaced starting from 1.02 (light) in 0.02 
decrements to smaller values (darker). 



phase that is expected to last for a few 100 ms and corresponds to a significant fraction 
of the overall signal, the u x flux may be almost a factor of 2 smaller than the v e flux, 
but it will be slightly hotter and less pinched [IB]. This corresponds to Case (a) above; 
it is evident from Fig. |3] and Table 0] that this implies that the Earth effect is very 
small. During the Kelvin-Helmholtz cooling phase the flux ratio is reversed with more 
v x being emitted than z/ e , but still with the same hierarchy of energies. This corresponds 
to Case (b); in this case the Earth effect could be about 6%. This time dependence may 
allow one to detect the Earth effect without a precise absolute detector calibration. 

In order to illustrate the time dependence of the Earth effect we show in Fig. [S]thc 
expected counting rate in IceCube for both the Livermore (left panels) and Garching 
(right panels) simulations. In the upper panels we show the expected counting rate 
with flavor oscillations in the SN mantle, but no Earth effect (solid lines), or with the 
asymptotic Earth effect (dashed lines) that obtains for a large Earth-crossing path. 
Naturally the differences are very small so that we show in the lower panels the ratio of 
these curves, i.e. the expected counting rate with/without Earth effect as a function of 
time for both Livermore and Garching. While for the Livermore simulation there is a 
large Earth effect even at early times, the change from early to late times in both cases 
is around 4-5%. Therefore, the most model-independent signature is a time variation 
of the Earth effect during the SN neutrino signal. 

In order to demonstrate the statistical significance of these effects we integrate 
the expected signal for both simulations separately for the accretion phase and the 
subsequent cooling phase; the results are shown in Table |3 For both simulations the 
Earth effect itself and its change with time is statistically highly significant. Based on 
the Livermore simulation, the Earth effect is much more pronounced than in Garching, 
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Figure 5. Earth effect in IceCube. The upper panels show the expected counting rate 
based on the Livermore (left) and Garching (right) models, including flavor oscillations. 
The solid line is without Earth effect (L = 0), the dashed line with asymptotic Earth 
effect (L ^ 0). The lower panels show the ratio between these curves, i.e. the ratio of 
counting rates with/without Earth effect. 



the latter involving more up-to-date input physics for neutrino transport. However, the 
difference between the Earth effect during accretion and cooling is not vastly different 
between the two simulations. Recalling that the absolute detector calibration may be 
very uncertain so that one has to rely on the temporal variation of the Earth effect, the 
difference between Livermore and Garching becomes much smaller. We expect that it is 
quite generic that the temporal change of the Earth effect is a few percent of the overall 
counting rate. 

Table 6. IceCube Cherenkov counts for the numerical SN models. 

Livermore Garching 





Accretion 


Cooling 


Accretion 


Cooling 


Integration time [s] 


0-0.500 


0.500-3 


0-0.250 


0.250-0.700 


SN Signal [Counts] 










No Earth Effect 


519,080 


818,043 


173,085 


407,715 


Asymptotic Earth Effect 


488,093 


751,137 


171,310 


390,252 


Difference 


30,987 


66,906 


1,775 


17,463 


Fractional Difference 


-5.97% 


-8.18% 


-1.03% 


-4.28% 


Background [Counts] 


720,000 


4,320,000 


360,000 


648,000 


•y/Background/Signal 


0.16% 


0.25% 


0.35% 


0.20% 
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3.3. Super- or Hyper- Kamiokande and IceCube 

One can measure the Earth effect in IceCube only in conjunction with another high- 
statistics detector. We do not attempt to simulate in detail the SN signal in this other 
detector but simply assume that it can be measured with a precision at least as good 
as in IceCube. One candidate is Super- Kamiokande, a water Cherenkov detector that 
would measure around 10 4 events from a galactic SN at a distance of 10 kpc. Therefore, 
the statistical precision for the total neutrino energy deposition in the water is around 
1% and thus worse than in IceCube. Even though Super-Kamiokande will measure a 
larger number of Cherenkov photons than IceCube, a single neutrino event will cause an 
entire Cherenkov ring to be measured, i.e. the photons are highly correlated. Therefore, 
in the estimated statistical y/~N fluctuation of the signal, the fluctuating number N is 
that of the detected neutrinos. If the future Hyper-Kamiokande is built, its fiducial 
volume would be about 30 times that of Super-Kamiokande. In this case the statistical 
signal precision exceeds that of IceCube for the equivalent observable. 

We denote the equivalent IceCube signal measured by Super- or Hyper-Kamiokande 
as N SK and the IceCube signal as iV Ic . If the distances traveled by the neutrinos before 
reaching these two detectors are different, the Earth effect on the neutrino spectra may 
be different, which will reflect in the ratio N^k/Niq. Of course, in the absence of the 
Earth effect this ratio equals unity by definition. 

The geographical position of IceCube with respect to Super- or Hyper-Kamiokande 
at a latitude of 36.4° is well-suited for the detection of the Earth effect through 
a combination of the signals. Using Fig. El we can already draw some qualitative 
conclusions about the ratio N SK /Ni C . Clearly, N SK /Ni C = 1 if neutrinos do not travel 
through the Earth before reaching either detector. If the distance traveled by neutrinos 
through the Earth is more that 3000 km for both detectors, the Earth effects on both 
iVsK and Nic are nearly equal and their ratio stays around unity. If the neutrinos come 
"from above" for SK and "from below" for IceCube, or vice versa, the Earth matter 
effect will shift this ratio from unity. 

In Fig. El we show contours of iVgx/iVic for the SN position in terms of the location 
on Earth where the SN is at the zenith. The map is an area preserving Hammer- Aitoff 
projection so that the sizes of different regions in the figure gives a realistic idea of 
the "good" and "bad" regions of the sky. In order to generate the contours we use 
the parameters of Case (b) in Table 01 so that the asymptotic suppression of the signal 
is about 6.5%. The sky falls into four distinct regions depending on the direction of 
the neutrinos relative to either detector as described in Table [3 When the neutrinos 
come from above for both detectors (Region D) there is no Earth effects. If they come 
from below in both (Region C), the Earth effect is large in both. Depending on the 
exact distance traveled through the Earth, the event ratio can be large, but generally 
fluctuates around 1. In the other regions where the neutrinos come from above for one 
detector and from below for the other (Regions A and B) the relative effect is large. 
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Table 7. Regions in Fig. for the Earth effect in IceCube and Super-Kamiokande. 

Region Sky fraction Neutrinos come from Nqk/Nic 

IceCube Super-K 

A 0.35 below above 1.070 

B 0.35 above below 0.935 

C 0.15 below below Fluctuations around 1 

D 0.15 above above 1 



4. Conclusions 

For assumptions about the flavor-dependent SN neutrino fluxes and spectra that agree 
with state-of-the-art studies, the Earth matter effect on neutrino oscillations shows 
up in the IceCube signal of a future galactic SN on the level of a few percent. If 
the IceCube signal can be compared with another high-statistics signal, notably in 
Super-Kamiokande or Hyper-Kamiokande, the Earth effect becomes clearly visible as 
a difference between the detectors. As one is looking for a signal modification in the 
range of a few percent, the absolute detector calibration may not be good enough in 
one or both of the instruments. However, for typical numerical SN simulations the 
effect is time dependent and most notably differs between the early accretion phase and 
the subsequent neutron star cooling phase. Therefore, one would have to search for a 
temporal variation of the relative detector signals of a few percent. The large number of 
optical modules in IceCube renders this task statistically possible. In fact depending on 
the differences in flavor-dependent fluxes, the statistical accuracy of Super-Kamiokande 
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may turn out to be the limiting factor. This limitation is not significant for Hyper- 
Kamiokande. 

The unique location of IceCube in Antarctica implies that for a large portion of the 
sky this detector sees the SN through the Earth when Super- and Hyper-Kamiokande 
sees it from above, or the other way round, i.e. the chances of a relative signal difference 
between the detectors are large. If both detectors were to see the SN from above there 
would be no Earth effect to detect. 

Assuming that the magnitude of the mixing angle #13 can be established to be large 
in the sense of sin 2 9 13 > 10~ 3 by a long-baseline experiment [321 HUH EHI > observing the 
Earth effect in SN anti-neutrinos implies the normal mass hierarchy. On the other hand, 
if sin 2 9i3 < has been established, the Earth effect is unavoidable. Not observing 
it would imply that the primary SN neutrino fluxes and spectra are more similar than 
indicated by state-of-the-art numerical simulations. 

If sin 2 # 13 >10~ 3 is known, and we do not observe the Earth effect, it still 
does not prove the inverted mass hierarchy. It could also mean that we do not 
properly understand the flavor-dependent source fluxes and spectra. Therefore, even 
if sin 2 #13 > 10~ 3 is known, our method only allows one to detect the normal mass 
hierarchy, it does not strictly allow one to exclude it. As far as neutrino parameters 
are concerned, only a positive detection of the Earth effect would count for much. Of 
course, a normal mass hierarchy and sin 2 #13 > 10~ 3 is certainly a plausible scenario so 
that expecting a positive identification of the Earth effect is not a far-fetched possibility. 

In summary, even though galactic SNe are rare, the anticipated longevity of 
IceCube and the long-term neutrino program at Super- or Hyper-Kamiokande imply 
that detecting the Earth effect in a SN neutrino burst is certainly a distinct possibility. 
This could identify the normal neutrino mass hierarchy, a daunting task at long-baseline 
experiments [321 ESI EH]- Given the difficulty of pinning down the mass hierarchy at long- 
baseline experiments, both IceCube and Super- or Hyper-Kamiokande should take all 
instrumental and experimental steps required to ensure the feasibility of a high-statistics 
simultaneous SN observation. 
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